Dispersion and Symmetry of Bound States in the Shastry-Sutherland Model 



O 

o 
o 

(N 



ON 



i 

c3 



I 

c 

o 
o 



> 

(N 
(N 

co 
in 
o 
o 
o 

a 

c 

o 
o 



X 
S3 



Christian Knetter, Alexander Biihler, Erwin Muller-Hartmann, and Gotz S. Uhrig 
Institut fur Theoretische Physik, Universitat zu Koln, D-50937 Koln, Germany 

(February 1, 2008) 

Bound states made from two triplet excitations on the Shastry-Sutherland (ShaSu) lattice are 
investigated. Based on the perturbative unitary transformation by flow equations quantitative 
properties like dispersions and qualitative properties like symmetries are determined. The high order 
results (up to ( J2/J1) 14 ) permit to fix the parameters of SrCu2(BOa)2 precisely: J\ — 6.16(fO)meV, 
x := J2/J1 = 0.603(3), J± = 1.3(2)meV. At the border of the magnetic Brillouin zone (MBZ) a 
general double degeneracy is derived. An unexpected instability in the triplet channel at x — 0.63 
indicates a first order transition towards a triplet condensate, related to classical helical order. 

PACS numbers: 75.40. Gb, 75.30.Kz, 75.50.Ee, 75.10.Jm 



Dedicated to Prof. F. Wegner on occasion of his 60 th birthday. 

H(J 1 ,J 2 ) = J 1 SiS J 

(ij) dimer 



Quantum antiferromagnets are at the center of re- 
search not only because of the high T c superconductors . 
Of particular interest are systems which do not have an 
ordered, Neel-type ground state. Their ground state is a 
spin liquid without long range spin order. Spin liquids 
are favored by low spin (S = 5 mostly), low coordina- 
tion number (Z G {2,3,4} =>• D e {1,2}), and strong 
geometric frustration. 

Dimcr solids are transparent cases of spin liquids. In 
D = 1, the generic example is the Majumdar-Ghosh 
model (!]] of which Shastry and Sutherland found a D = 2 
generalization (ShaSu model) In both cases frustra- 
tion is essential. Each spin is coupled to pairs of spins 
(dimers). If these pairs form singlets the couplings be- 
tween dimers is without effect and the singlet-on-dimers 
product state is always an eigen state and for certain 
parameters the ground state lU^H . The systems are 
gapped. The elementary excited states are dressed 5 — \ 
(D = 1) H or S = 1 (D = 2) entities. They inter- 
act strongly and form bound and antibound states in 
various spin channels. Due to its recent realization in 
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FIG. 1. Unit cell of the Shastry-Sutherland (ShaSu) lattice 
as realized in SrCu2(B03)2, (a,b) are unit vectors; (a',b') 
and the coefficients 7 are used to analyse Raman scattering. 



SrCu 2 (B0 3 )2 §0 the ShaSu model (Fig. fjj) is presently 
attracting enormous interest. The Hamiltonian reads 



In this Letter we start from the dimer phase ||. We 
focus on bound states formed from pairs of the elemen- 
tary triplets and their symmetries, degeneracies, and dis- 
persion. The perturbative unitary transformation j^] 
based on flow equations || enables us to link smoothly 
and uniquely H(Ji, J 2 ) at x := J2/J1 7^ to an effec- 
tive H c s conserving the number of triplets on dimers 
[H e ff, H(Ji, 0)] = 0. This permits a clear distinction be- 
tween the ground state (without triplets), the 1-triplet 
sector, the 2-triplet sector etc.. 

In terms of H e g the dynamics of one triplet is hopping 
th-.i (ty.i) starting from a horizontal (vertical) dimer by i x 
dimers right and i y dimers up. Nothing else is possible 
due to triplet number conservation. The elements t are 
computed in order 15 |3 10 T^. 

The dynamics of two triplets at large distances is gov- 
erned by 1-triplet hopping. At smaller distances a 2- 
particle interaction occurs additionally given by Wh t d t i,d' 
(W^;d;i,d') starting with one triplet on a horizontal (ver- 
tical) dimer and another at distance d. The action of 
H c s is to shift the triplets to i and to i + d'. Nothing 
else is possible due to triplet number conservation. Since 
the total spin is conserved (S S {0, 1,2}) the distances 
are restricted to d, d' > 0, i.e. d x > or d x — A d y > 0, 
because the exchange parity is fixed. 

The action of H c fi yields the combined effect of hopping 
and interaction denoted by Ad;i,d'- The true 2-triplet 
interaction is easily found by subtracting the 1-triplet 
hopping §00 

Wd-0,d' — ^d;0,d' — £d'-d — ^d'.d^O (2a) 
W^d-d'.d' = ^4d;d-d',d' — ^d-d' — ^d'.d^O (2b) 

Wd ; -d',d' = A d; _ d /,d' - t-d-d' (2c) 
W d ; d ,d' = A d; d,d> - t d +d' (2d) 

(distinction h/v omitted for clarity). Otherwise A and W 
are identical. The coefficients W for S £ {0, 1, 2} yield 
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the complete 2-particle dynamics. We compute W up to 
x 12 . the coefficients for the lowest-lying states even up to 

During the virtual processes Jl2| the triplet number is 
changed. Due to the frustration of the ShaSu lattice H 
in (Q) changes the number of triplets on the dimers at 
most by one j^,|ll|[ll|] . An excitation or a de-excitation 
on a horizontal (vertical) dimer is possible iff at least one 
of the vertical (horizontal) dimers on the left and right 
(above and below) are excited. This restriction implies 
that one triplet hops only in x 6 (cf. Figs, in [ [l3||i"T| | ) . 

Motion of two triplets together is much less restricted 
(cf. Fig. H). Matrix elements occur in x 2 as first ob- 
served for total spin S = 2 But the dispersion of 
bound states starts only in x 3 (contrary to x 4 claimed 
in Ref. fl5| ). Two adjacent triplets interact linearly in x 
(-x for S = 0, -x/2 for S = 1, x/2 for S = 2). The 
energy of the initial and final state in each row in Fig. || 
differ by 0(x). Hence both rows must be combined mak- 
ing it an (x 2 ) 2 /x = x 3 process eventually. This applies 
to the 8 (anti)bound states derived from two triplets on 
(next) nearest neighbor dimers. The dispersion of any 
other state sets in at higher order. 



x 2 process ! 



x 2 process ! 



FIG. 2. Leading processes of correlated 2-triplet hopping. 
Dark dots are triplets, bars are dimers. 



We use the following basis for the 2-triplet states 




where k is the conserved total momentum in the magnetic 
Brillouin zone (MBZ) applying due to the two sublattices; 
a e {0,1}, Q := (tt, 7t), N is the number of dimers, 
r, r + d) denotes the state with triplets at r and at r + d. 
The distance d is restricted d > 0, i.e. d x > or d x — 
A d y > 0. The matrix elements of H e g in the basis (j^) 
are real due to translation and inversion symmetry. 

Before the quantitative analysis of the bound states 
a qualitative aspect, a general double degeneracy at the 
border of the MBZ, shall be derived. To see this consider 
the combination of a shift by the dimer-dimer spacing 
along a' (S), a reflection about a (R), and the inversion 
r — > r (I) (cf. Fig. |l|). The combinations SR and I are 



symmetries of the Hamiltonian. For k x 
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|k,d,£7> -> e^ + "^ +ti » + ' T+PS )|k,-( J ),1 -a) , (4) 

where P <E {0, 1} being unity iff — (d y ,d x ) < so that 
the triplets must be swapped to pass from — (d y ,d x ) to 
(d y ,d x ). It is crucial that SRI links \k,(d x ,dy),0) and 



|k, — (dy, d x ),l) like a 2D rotation 
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up to a pref- 



k, d, a) := N- 1 ' 2 V e Mk+ CT Q)(r+d/2)) ^ r + d) ^ (3) ^ 



actor. Hence its eigen vectors are complex with linearly 
independent real and imaginary part and so are the si- 
multaneous eigen vectors of SRI and H Q g . Because H c s is 
real the real and the imaginary part constitute in fact lin- 
early independent eigen vectors to the same eigen value. 
The same double degeneracy is concluded for the other 
parts of the MBZ border by S and 90° rotation (D). It is 
also valid in the 1-triplet sector pdj |. 

The double degeneracy at the MBZ border is interest- 
ing for analysing experiment, too. Degeneracy reduces 
the large number of energetically close states helping to 
resolve different bound states. 

Since 1-triplet hopping is of higher order than in- 
teraction an analytic expansion for the energies of the 
bound states is possible. At finite order in x only con- 
figurations contribute where the two triplets are not 
too far away from each other. Of course, higher or- 
ders imply larger, but still finite distances. In par- 
ticular, the energies of the four states which evolve 
from neighboring triplets can be computed very well 
since their interaction is linear. Investigating the ma- 
trix elements shows that it is sufficient to study the dis- 
tances d e {(0,1),(1,0),(1,±1)} for order 5. To x 14 
only d £ {(1, ±2), (2, ±1), (0, 2), (2, 0), (2, ±2)} must be 
added. So, for given total momentum only a finite 8x8 
or 24 x 24 matrix has to be analysed. For illustration 
consider the elements ^4(o,i) ; i.(2,i) (the Fourier transform 
of i yields the momentum dependence.) connecting (0, 1) 
and (2, 1) which is 0(x 4 ). By second order perturbation 
one sees that the resulting energy shift is (x 4 ) 2 /x = x 7 



of the MBZ border) definition ([|) implies for the total 
combination SRI the mapping 



Furthermore, the elements connecting shorter dis- 
tances to longer distances and the elements among longer 
distances do not need to be known to very high order. 
Consider again the process (0,1) <-> (2,1). In order x 7 
the element -A(o.i) ; i,(2.i) must be known only in x 4 and 
^(2.i) ; i,(2.i) on ly i n x ; in order x 9 the element A(n,iy,U2,i) 
must be known only up to x 6 and ^4(2.i) ; i,(2.i) only in x 3 
and so on. 

We have analysed the dispersions in x 5 of the four 
states bound linearly in x in the MBZ. Fukumoto's re- 
sults are mostly confirmed Jl6|,|l7]]. At particular points 
of high point group symmetry ((0, 0),(0, n),(ir/2, n/2)) 
the Hamiltonian splits into six blocks corresponding to 
different representations of the square point group 4mm. 
At these points the analysis up to x 14 is carried out Jl7| . 
The symmetries are classified according to the irreducible 
representations (four ID, one 2D) of the point group 
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4mm IM1), T 2 (x 2 - y 2 ), T 3 (xy), r 4 (^(x 2 - y 2 )), T 5 (x, y) 
where simple polynomials are given in brackets to show 
the transformation behavior. 

The extrapolated energies are depicted in Figs. || (5 = 
0) and [| (5 = 1) as functions of x. For those ener- 
gies which stay separated from the 2-particle continuum 
Dlog-Pade approximants are used successfully The 
results are stable under changes of the polynomial de- 
grees. The energies close to the continuum (here simply 
twice the gap A between the ground state and the ele- 
mentary triplet at k = 0) are given with less reliability 
by the truncated series or by a non-defective Dlog-Pade 
approximant. 




0.3 0.4 0.5 0.6 0.7 

x=J 2 /J 1 

FIG. 3. Energy of the lowest-lying S — states. Curves 
refer to k = except the dashed-dotted one. The dotted 
curve displays the continuum at 2A. Inset: 1-triplet disper- 
sion. Theory at x = 0.603, Ji = 6.16meV; data from Ref. 
experimental errors at least as large as indicated. 



In Figs. ||, H the modes are sorted in energetically as- 
cending order for small values of x: solid, long dashed, 
and short dashed curves. The Y§ modes are naturally 
degenerate. The double degeneracy for k = (0,7r) does 
not result from the point group but originates from the 
complex conjugation as explained above. The dashed- 
dotted curve at (0, n) has to be compared to the solid 
and the long-dashed curve to assess the dispersion of 
these two modes from to (0,7r). While for S = 
this dispersion always has the expected behavior with 
to(0) < w((0, 7r)) the energies for S = 1 are reversed for 
small values of x (cf. [jl6|). Only above x m 0.55 the 
relation w(0) < w((0,7r)) holds for 5=1. 

We do not agree with Ref. (T^| that the two lowest 
states are of s-wave type since this would imply that 
they are Ti. Instead the 5 = 0, 1 states are odd under 
reflection about a' (Rl) or about b' (R2) (cf. Fig. 0). For 
5 = 0, the lowest state is even under SD and the second 
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x=J 2 /J 1 

FIG. 4. As in Fig. ^ for 5 = 1. Inset: Magnetic suscepti- 
bility. Theory (dashed) with directional rms average g = 2.13 
po[ , x, Ji as in Fig. |^; experiment (solid) on powder |H|. 



lowest is odd. For 5 = 1, it is vice- versa. The T$ states 
can be viewed as being of p-wave type. 

For 5 = 0, the lowest mode vanishes at the same x 
as does the gap A. So no additional instability occurs 
for 5 = 0. This provides evidence against a competing 
singlet phase as presumed in Ref. pT| . There is, how- 
ever, a salient instability for 5 = 1 (Fig. |) at i = 0.63. 
This comes as a surprise since one expects in antiferro- 
magnets binding effects to be largest for 5 = 0. The 
singularity at x = 0.630(5) is very stable occuring in all 
non-defective Dlog-Pade approximants. We take the van- 
ishing of a bound 2-triplet state with 5 = 1 as evidence 
for strong attraction between triplets which are neither 
parallel (5 = 2) nor antiparallel (5 = 0). The attraction 
points towards a first order transition into a condensate 
of triplets occuring at much lower x than thought previ- 
ously |g[ . The angle being neither zero nor it corroborates 
a helical order if one adopts a classical view |^,|| . In the 
light of the instability for 5 = 1 we interprete the findings 
by Koga and Kawakami [2l[] as indication of the same in- 
stability. It is so far not excluded that the bulk triplet 
condensate (the helical phase) is a singlet or that it can 
be linked to a singlet p3[ . 

Next we determine x and Ji for SrCu2 (BOa)2. The 
steep decrease of the bound 5 = 1 state enables us to 
fix x very precisely. Based on ESR (20), FIR 24| as 
well as INS [jl9| we assume A = 2.98meV and w|s=i = 
4.7meV leading to x = 0.603(3) and J x = 6.16(10)meV. 
The 1-triplet dispersion is in excellent agreement with 
experiment (cf. inset in Fig. |^ and Ref. |ll]]). Raman 
scattering |l5| provides further strong support because 
the energy of the T3 singlet matches 30cm _1 perfectly 
The T 4 singlet at 25cm -1 [p5| is forbidden by symmetry 
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since the Raman operator is effectively T3. 

In leading order t/U the Raman operator R = 
^7ijSjSj couples the same spins as the Hamiltonian. 
But only the antisymmetric (71 = —72) part of R on 
the dashed bonds (cf. Fig. |l|) creates excitations from 
the ground state. By geometry we have 7" = 72 and 
7 2 ' = 71 so that the effective i? c ff is odd under Rl and 
R2. But the projection of the vector potential A (A||E) 
on the bonds under study in the polarisation (ab) implies 
7i = 72 = — 7i = —72 = 7Cos(2a) (7 microscopic con- 
stant, a angle in Fig. |l|), i.e. an even component imply- 
ing R c ff = 0. On the contrary, polarisation (a'b') yields 
71 = —72 = — 71 = 7 2 = 7 sin(2a) so that i? G ff ^ 0. This 
finding agrees nicely with experiment where on T — > the 
intensities almost vanish for (ab) but grow for (a'b') [fj"5| . 
Additionally, 71 = —7^ and 72 = — 7 2 imply odd parity 
under SD so that R c s is indeed T3, not T4. Calculat- 
ing the next r 3 , S = bound state (less systematically) 
yields 45cm~ 1 in good agreement with the experimental 
46cm~ 1 line, too. 

We conclude that the 2D model ([!]) explains the low- 
lying excitations of SrCu2 (603)2 perfectly. Thermody- 
namic quantities like the susceptibility x{T) require the 
inclusion of the interplain coupling J± which is fully 
frustrated not changing the dimer spins [ p6| . We have 
employed a Dlog-Pade approximant for the 1/T expan- 
sion of the 2D X2D JR} complemented by the condition 
A = 2.98meV. This ansatz works fine for T > 35K. The 
3D X3B is computed from \2D on the mean-field level as 
X3D = X2D + 4Jj_ ■ The inset in Fig. shows that theory 
(J_L = 1.3(2)meV) and experiment M agree without flaw 
above 40K. Our value for J± is significantly higher than 
the one in Ref . [^6| due to different values of x and J\ . 

The above comprehensive analysis of bound states is 
a fine example of the efficiency of perturbation by uni- 
tary transformation [0 based on flow equations. This 
clear concept allows to distinguish uniquely sectors with 
different particle numbers and other different quantum 
numbers like the total spin. Here the concept was put 
to use to analyse the Shastry-Sutherland lattice as real- 
ized in SrCu2(B03)2- To our knowledge it is the first 
quantitative description of 2-particle bound states in 2D. 

The symmetries of experimentally relevant states were 
determined. The reliability of the high order results 
allows to fix the experimental coupling constants very 
precisely (Ji = 6.16(10)mev, J 2 /J\ = 0.603(3), J± = 
1.3(2)meV). Thereby, different experiments (ESR, FIR, 
INS, Raman, x(T)) are explained consistently We sug- 
gest to exploit the double degeneracy derived here to re- 
solve different bound states at the border of the MBZ. 

An unexpected instability for the S = 1 2-triplet bound 
state is found at x « 0.63 indicating a transition to a 
triplet condensate probably related to the helical phase 
found previously p2|,p[. We conjecture that this tran- 
sition is first order occuring at lower x than assumed 
so far. The vicinity of SrCu2(B03)2 to this transition 



suggests to attempt a closer experimental analysis. Pres- 
sure and/or substitution will certainly influence the ratio 
J2/J1. Thereby one may hope to scan through the tran- 
sition and to examine the phase beyond. 

The authors like to thank H. Kageyama and P. Lem- 
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and discussion. The work is supported by the DFG in 
SFB 341 and in SP 1073. 
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